Complete positivity and contextuality of quantum dynamics 
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Positivity or the stronger notion of complete positivity, and contextuality are central properties of quantum 
dynamics. In this work, we demonstrate that a physical unitary-universe dilation model could be employed to 
characterize the completely positive map, regardless of the initial correlation condition. Particularly, the prob- 
lem of initial correlation can be resolved by a swap operation. Furthermore, we discuss the physical essence of 
completely positive map and highlights its limitations. Then we develop the quantum measurement-chain for- 
mula beyond the framework of completely positive map in order to describe much broader quantum dynamics, 
and therein the property of contextuality could be captured via measurement transfer matrix. 



I. INTRODUCTION 

Quantum dynamics, which lies in Hilbert space, is ex- 
plicitly distinct from classical mechanics and stochastic dy- 
namics, although there are some similarities and connections 
among them. As the properties of quantum dynamics, su- 
perposition and uncertainty have been explored greatly, while 
nowadays along with the development of quantum control and 
computation |l|,|2|] etc, the notions of complete positivity, en- 
tanglement, nonlocality, and also contextuality are of funda- 
mental concern and more interest to researchers. 

The generic frameworks to describe quantum dynamics in- 
clude theories of quantum channel dH 0] and quantum tra- 
jectory [5], in this work we focus mainly on the former one. 
As quantum state, namely density matrix, could be treated 
as one positive semi-definite operator, quantum dynamics is 
characterized as a map acting on the state. A formal theory 
based on positive and completely positive (CP) maps provides 
a well-defined mathematical framework [3, 4], from which the 
microscopic Markovian dynamics can be derived associated 
with a semigroup structure iH-IH]. The CP map, which is also 
known as a quantum channel, corresponds to the operator-sum 
representation (OSR) [9]. However, the limitation is that CP 
map is not enough to describe all quantum dynamics, since 
CP map is only a subset of positive map. Recently the lim- 
itations of CP map cause much attention, namely, the initial 
correlation issue. It is proved that a map with initial classical- 
quantum state, which has null quantum discord, can be written 
as a CP map II 1 0L 1 1 111 . We will show that, via a swap operation 
trick, the open-system dynamics can always be written as a 
CP map when there is any initial correlation. 

We analyze the physical properties of CP map employing 
a "unitary-universe dilation model". We then show that the 
dilation model has limitations when applied to describe all 
quantum dynamics, for which we rather explore a different 
framework for the rest of the work. It is proved that, beyond 
the dilation model, the OSR form with state-dependent Kraus 



operators could be constructed for a positive, non-CP map 
between any two quantum states IIX2h . The state-dependent 
OSR form could be understood as state-dependent quantum 
measurement with memories and feedbacks. One construc- 
tive and operational quantum measurement-chain formula is 
developed in this paper to describe quantum dynamics be- 
yond the dilation model. The quantum measurement trans- 
fer matrix, which is the analogue of transfer matrix in clas- 
sical stochastic dynamics, captures the contextuality of quan- 
tum dynamics. Contextuality, originating from the Kochen- 
Specker theorem lfl3ll . can be described via inequalities lfl4tl 
or measures of contextuality [15]. A common example of 
contextuality of quantum observable is that for observable A 
commuting with B and C, which do not commute with each 
other, the measurable values of A can be different when it is 
measured together with B or C. Technically, the reason is 
that one operator can have different decompositions, i.e. fine- 
grained structures, which can lead to physical effects in prac- 
tice. In this work we take a generalized notion of contextual- 
ity [16], without focusing on the exact distinctions of different 
notions. Briefly, we say there exists contextuality of one op- 
erator (density matrix, observable, POVM (positive operator- 
valued measurement) etc), if the fine-grained structure of the 
operator causes measurable physical effects. 
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In this work, we study issues relating to properties of com- 
plete positivity and contextuality, by employing the unitary- 
universe model and quantum measurement-chain formula, re- 
spectively. The paper mainly contains two parts. In section HT1 
we study the mathematical framework for CP map, and par- 
ticularly analyze the initial correlation problem, which can be 
solved via a swap trick in the unitary-universe dilation model. 
We define a physically positive map to address the limitations 
of CP map. Also, the limitations of Bloch vector dynamics 
in Heisenberg picture are discussed. In section [HI] we study 
the quantum measurement-chain formula, and properties of 
measurement and weak-measurement transfer matrices. The 
contextual property of quantum dynamics is demonstrated. At 
last we conclude briefly in section HVl 
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FIG. 1 : Schematic diagram for the unitary-universe dilation model. 
U and W represent unitary operations. The meters represent trace 
operations. We introduce in the notation for convenience, that if two 
systems S and E are correlated, we denote S ooE, and if two systems 
5 and E interact (do not interact), we use 5 «-» E (S m E). 



n. COMPLETELY POSITIVE MAP 

In this section, we briefly review the standard framework 
of CP map, and then establish the unitary-universe dilation 
model to interpret the physics. For the initial correlation prob- 
lem, we employ the swap operation and discuss one example 
to demonstrate its viability. Then we discuss the properties 
of physically positive map and hermitian map, and highlight 
their limitations. 

According to operator theory, quantum state p is treated as 
a trace class, positive semidefinite operator acting on Hilbert 
space H, and p€f (H). A map Q> : T (H) — > T (H) is positive 
iff it maps a state p into a state <J>(p). A map <t> is completely 
positive (CP) iff O® 1,„ is positive for any m e Z + . Stinespring 
dilation theorem and Choi's theorem llH |H 03] state that a CP 
map takes the form 



(1) 



!=1 



and the dilation is said to be minimal for N = n , where n is 
the dimension of the Hilbert space H. If Yn=\ K ■ K, = 1, the 
map is trace-preserving (TP). Note that a CP map <I> induces 
the CP map <I> ® 1,„, i.e. if map <I> is CP, then for all m, the 
positive map <& ® 1,„ is also CP, which follows from the tensor 
product property of CP map. 



A. Unitary-universe dilation model 

Kraus studied the physical essence of CP map, and es- 
tablished the operator-sum representation (OSR) [9], which 
forms a physical dilation model for CP map. The model we 
study below is termed as unitary-universe dilation model (u- 
u), in which a system S is distinct from the rest of the universe 
E, and the universe evolves unitarily. In the usual case, sup- 
pose the initial universe state is the product of system state 
ps and environment state p E , and the global evolution is uni- 
tary U, then the final state of the system, denoted as p' s , after 
tracing out environment is 



p's = ®(Ps) = tt E (U Ps ®p E U r ) 
= J^W>Ps<eiE/V>, 



where p £ = 2; M e ')( e '\> and the Kraus operator will take the 
form K,j = yfAi{e^\U\e'}. The minimal dilation is attained 
when the Kraus rank is n 2 for n dimensional system. 

The OSR form has the properties as follows: i) the Kraus 
operators are independent of system state; ii) the map cD ® 1,„ 
acting on psm could be written in OSR form, with M as arbi- 
trary ancilla acted on by 1,„ with dimension m; iii) there are 
unitary freedom for the Kraus operators. All the properties are 
consistent with the requirement of the dilation theorem, which 
is that there cannot be initial correlation between system and 
environment, since the state ps m needs to be arbitrary. In other 
words, the requirement of CP is in conflict with the initial cor- 
relation between system and environment. This point agrees 
well with Pechukas's argument 111 80 (on 3ed page) that it is 
not always possible to prepare a global state Psem, where sys- 
tem S and ancilla M are correlated while the reduced state of 
system and environment pse is as arbitrary as we choose. 

However, in physical reality, initial correlation between 
system and environment is possible. This is the main origin of 
the CP problem. The initial condition for u-u model is S ooE, 
and S E, as shown in the left panel of Fig. Q] The dimen- 
sions of S and E are arbitrary. 

There have been some results dealing with the initial corre- 
lation problem ITol [Til I20I - I23I1 . It is proved that a map with 
initial vanishing quantum discord type correlation can be writ- 
ten as CP map fiUl . However, the map will depend on the 
system state, and the system state can only be in a restricted 
class, thus conflicting with the requirement of CP. A positive 
map with initial correlation can be written in afRne map form 
Il22ll . yet it is not CP anymore. One approach is developed em- 
ploying one so-called assignment map Izjl , which needs to 
be linear, consistent, and positive at the same time. However, 
it is shown that the assignment map can only be linear at the 
expense of giving up positivity or the consistency. Despite all 
the difficulties, it turns out by a simple swap trick studied be- 
low, which also has been used in other studies lfll l24l - l27ll . the 
u-u model corresponds to a CP map on the system regardless 
of the initial condition. 

Suppose the state of system S acts on n dimensional Hilbert 
space H, and there exists a copy of this space which then 
forms a part of the environment Hilbert space, the correspond- 
ing system is denoted as ancillary-environment A. Let the 
state of S + E be p SE , the state of A be p A = ti E {p SE ), and 
the swap operator Vsa act between the system and ancillary- 
environment before the evolution of the universe, and thus 



Pse ®Pa = (V S a ® Was ® cr £A (V ® 1 £ ) 



(3) 



(2) 



where as = p A = tr £ (p S£ ), ct E a = Pse, and V^Vsa = 1- As 
the result, the dynamics with initial condition p$E ® Pa and 
evolution Use ® 1a is equivalent to that with initial condition 
os ®cr EA (i.e. m E (p SE )®p SE ) and evolution (U S e®1a){Vsa® 
1 E ) = W. The swapped u-u model is depicted as the right 
panel in Fig. Q] The final system state p' s could be expressed 
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as 

p' s = tv EA (Wo- s ®o- EA W v ) 

= Y j P^ea(Wcts ® |<A,)<(A/|W + ) 

i 

= E^4' (4) 

hi 

where Kraus operator Kjj = ~\[pi(if/j\W\i[/i), and eigenvalue 
decomposition ct E a = E; Pilftdi^il- The form (0|i defines a 
CP map which maps the system state cr s = ti'EipsE) into p' s 
by the set of Kraus operators {Kjj}. Note in the above the 
ancillary-environment A is different from the ancilla M in the 
discussion of CP map. The choice of the initial state of the 
ancillary-environment A is unique since for the case with ini- 
tial correlation, given the state p$E, the only reasonable initial 
system state is ti'EipsE), which is just encoded in the state of 
the ancillary-environment. 

As to the dimension issue of the swapped u-u model, it is 
clear to observe that if the total dimension of the Hilbert space 
for E+A is n 2 , the model becomes rightly the minimal dilation 
model for CP map. In this case, the dimension of the Hilbert 
spaces for E and A is both n. In other words, we find that the 
swapped u-u model is indeed a particular kind of model for 
dilation, with bipartite environment and a special form of uni- 
tary evolution W. This restriction of dilation results from the 
initial correlation, since the formation of the initial correlation 
itself already reveals some information of the dynamics. 

We conclude that a CP map requires dilated unitary evolu- 
tion and null initial correlation, and can always be character- 
ized by the unitary-universe dilation model, regardless of the 
physical initial correlation condition. As the result, we extend 
the validity of CP map to the case with initial correlations pro- 
vided the whole universe in the model evolves unitarily. 



B. Convert non-completely positive map to completely positive 

map 

We analyze the form of swap trick in this subsection. If the 
system is a single qubit, the required swap operation is sim- 
ply the two-qubit swap gate SWAP = | 2f =0 o", ® cr,;. For a 
qutrit system, it is the two-qutrit ternary swap gate 112611 . For 
many -body open-system with initial correlation, the dynamics 
is more complicated generally, since there could be multipar- 
tite correlations and £-local interactions. A form for general 
qudit swap gate and the algorithm to design it has been con- 
structed 02711 . For practical applications, the circuit design and 
implementation of the swap operation are feasible. 

Next, we discuss one simple yet precise example to show 
that how to convert non-CP map to CP map. Suppose a sys- 
tem qubit initially forms the maximally entangled state with 

the environment qubit as = -^(|00) + |11», where the 
states |0, 1} can be polarization, electron, or phonon states etc 
in practice. The simple evolution on the two-qubit "universe" 
is assumed to be U = e^' 6X ® x , with 8 = gt, g as coupling 
constant, t as time. The dynamics of the single-qubit system 



cannot be written as CP channel generally due to initial en- 
tanglement. With the method of swap trick above, one ancil- 
lary qubit is added to the universe with maximally mixed state 
Pa = jdiag(l, 1). According to formula ©, and using the 
Bell's basis to cany out the trace, the CP map on the system 
qubit with state cr s = |diag(l, 1) will be represented by four 
Kraus operators corresponding to a minimal dilation, 

K x = -e- m l, K 2 = -e~ m X, 
2 2 

K,= l -e- ie Y, K A = l -e- ie Z, (5) 

which is a single-qubit depolarization channel. Another way 
to demonstrate the CP map is to show the positivity of the 
Choi's matrix Jj]. According to channel-state duality, a chan- 
nel <D on a qudit is CP iff the Choi's matrix C = <i(<D<g>l)(|£X£|) 
is positive, with maximally entangled state |£) = I")/ V^- 
It is direct to check that the Choi's matrix for the depolariza- 
tion channel in this example is the identity operator, which is 
clearly positive. 



C. Limitation of the unitary-universe model 

The only physical requirement of quantum dynamics is pos- 
itivity, and so is the tensor product of quantum dynamics (in 
Schrodinger picture). Yet, there exists mismatch between the 
positive map and physically realizable processes. The domain 
of a positive map is the whole set of state space T(H). While 
some physical operations, e.g., partial transpose, can only be 
realized on a restricted class of states, although they are not 
positive. All positive map is supposed to be realized physi- 
cally. However, the problem is whether all physical dynamics 
is positive or not. Let us define a new class of map, the phys- 
ically positive (PP) map. Observe that a physical dynamics 
may only live in a subspace of the whole Hilbert space and re- 
alize a restricted class of state. We can introduce the class of 
"non-universal positive" map, which acts on state p in a sub- 
space of T(H). Then the PP map contains the positive map 
acting on all state in T(H) and the" non-universal positive" 
map on a restricted class of state. As the result, we have the 
set relation 

H D PP dPd CP D M, 

where we have included the set of hermitian map H and 
Markovian dynamics M for completeness. Consequently, we 
can describe partial transpose as a PP map acting on the class 
of separable state plus states which also stay positive under 
partial transpose. 

Now let us analyze the limitations of CP map and the u- 
u model. In the u-u model, the unitary evolution U of the 
universe can be reduced to a map on system S as <t>, or a cor- 
responding map on environment E as fl>'. The map <t> which 
acts on system is CP iff O ® l m is positive for any m and an- 
cilla M. However, the physical problem is that the ancilla M 
needs to exist, which cannot be a part from the environment E, 
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FIG. 2: Schematic diagram for the quantum measurement-chain 
model. The quantum system is prepared in p , the measurements 
M„ are performed on the system resulting in states p n , n = 1, 2, • • • , 
the quantum transfer matrix is T. The arrow represents the direction 
of evolution. 

since the unitary operator U acts containing the whole envi- 
ronment. Thus the ancilla M can only be outside of this origi- 
nal universe of S + E. This point indicates that the essence of 
CP map is not whether there is initial correlation or not, but 
instead, it is that the universe needs to evolve unitarily. For 
the case with initial correlation, we extend the universe with 
the ancillary-environment A, and the new universe S + E + A 
still evolves unitarily such that ancilla M is outside. As the 
result, for quantum dynamics beyond CP, the universe has to 
evolve non-unitarily. Note that if the evolution of the universe 
is a CP channel, then it can be further dilated to a unitary 
evolution resulting a CP map on system again. For the non- 
unitary universe model, the map on the universe is restricted 
as <b v e PP \ CP, resulting in another map O s e PP \ CP 
on the system. However, there will not be a precise universal 
formula for PP map. 

III. QUANTUM MEASUREMENT-CHAIN FORMULA 



and T. Thus, X„ = T"Xq. Classical Markovian chain satisfies 
Chapman-Kolmogorov equation 

p(x 3 , f 3 |*i,fi) = J dx 2 p(x2, t 3 \x 2 , t2)p(x 2 , t 2 \xuh)- (6) 

This equation has an alternative form, pix^J^) = 
j dx 2 dxip(x 3 , t 3 \x2, ti)p(x2, f2lx1.fi). 

Quantum dynamics is expressed as matrix equations (for fi- 
nite dimensional case), and thus looks similar with stochastic 
dynamics. However, there are differences between quantum 
dynamics and classical stochastic dynamics. First, for the no- 
tion of probability, classically it originates from the ignorance 
of the detailed underlying causal dynamics, which means that, 
the uncertainty will be eliminated if the details are available. 
While probability in quantum dynamics comes from measure- 
ment and the related decoherence of superposition or entan- 
glement. Furthermore, based on OSR form, quantum mea- 
surement is actually not random, that is, in a measurement 
\Mj} (2, M. Mi < 1), each operation M,- results in exactly one 
corresponding result. The other difference is that quantum dy- 
namics contains coherence, and it is contextual, i.e., different 
basis can exist in reality leading to various decompositions; 
on the contrary, classically only one basis exists and thus 
the dynamics is decohered or de-contexted, for instance, the 
state X and transfer matrix T have a common basis, namely, 
{0, 1, 2, ... }. In appendix lAl we show that the channel matrix 
% for a dynamics and the quantum state can be expressed in 
the same operator basis; however, quantum dynamics is one 
affine map instead of a linear map, and the channel matrix % 
is not a probability matrix. 



In this section, we investigate a theoretical framework dif- 
ferent from the CP map. For comparison, the Bloch coher- 
ence vector and affine map form is discussed in appendix lAl 
where its limitations are analyzed. There are basic differences 
among the various frameworks to be highlighted. The affine 
map form takes the quantum state as a vector, which can be 
visualized in the configuration space. However, the coherence 
vector is generally multi-dimensional, and the dynamics of 
many observables need to be determined in Heisenberg pic- 
ture. While in CP map form, the quantum state is treated as 
a whole, i.e. one operator, so the detailed structure (e.g. dif- 
ferent decompositions) of the state cannot be revealed. As we 
know, although quantum state is a positive operator, the same 
density operator could be physically prepared differently. In 
other words, the fine-grained structure of quantum state has 
physical consequences, which manifests the contextual prop- 
erty 11611 . By analogue with classical probability theory, below 
we develop a different technique to describe quantum dynam- 
ics, and the property of contextuality can be directly captured. 

Let us recall some basic classical probability theory $2^. 
The classical system is represented as a random variable X, 
which at each particular time t can take values within a set. 
A single-step (one-shot) stochastic dynamics is described by 
transfer matrix T. A stochastic dynamics is Markovian iff T 
is independent of variable X, and X„ only depends on X„_j 



A. Quantum measurement transfer matrix 

A quantum measurement-chain is a chain of measurements 
acting on the state p of a quantum system. Denote the chain 
with« measurements as M\^„ - M„oM„-\ °- • -oM\, which 
is depicted in Fig. [2] A quantum dynamics can be modeled 
as a measurement-chain, since any two quantum states can 
be connected via the OSR form 01211 . The measurements do 
not need to be of the same rank, i.e. the number of operations. 
The effect of the measurement-chain is just the composition of 
their corresponding channels. Although state p„ after the n-th 
measurement only depends on the measurement M„ and p„_i, 
the measurement M„ can contain the information of the states 
of former steps and the initial state. Thus the measurement- 
chain is generally non-Markovian. A standard Markovian dy- 
namics is described by Lindblad equation, which can be de- 
rived from the CP map |6j-|8|] with state-independent and time- 
independent Kraus operators. We say a quantum dynamics is 
non-Markovian if it can not be described by Lindblad equa- 
tion and equivalently other forms such as Redfield equation. A 
measurement-chain becomes Markovian if each measurement 
in the chain is independent of state and occupies infinitesimal 
time; then the generator of the channel is the sum of those for 
the building channels via Lie-Trotter decomposition Il28ll . In 
addition, there are different quantum Markovian-chain theo- 
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ries in literatures, e.g. quantum walk theory [j29j, and consis- 
tent history theory OOll . yet there is no measurement-chain. 

Let us construct a measurement-chain with two measure- 
ments M = {M/} and N = {Nj} for the sake of convenience, 
and further generalization is straightforward. Note the mea- 
surements could depend on system state generally. Suppose 
there is no evolution except during the measurement pro- 
cesses. A quantum state prepared in po, after the first mea- 
surement M, becomes 



Pi 



= Y,M;p M' 



(7) 



which can also be written as p\ = p/tr,-, with pi = 
tr(MJMipo), YiiPi - 1- After the second measurement AT, 
the state becomes 



P2 



J>pi4 



(8) 



with qj = tr(N'jNjPi), YijQj = !• The probability for the 
second measurement can be decomposed as 

qj = tr(.Y ; .\>i.) 

= J^HNjNjMiPoMj) 



ti(M N NjM iPQ ) 
X \ HMjMiPo) 



(9) 



rrijj forms the quantum measurement transfer matrix T, which 
is the analogue of classical transfer matrix T. The quantum 
states after each measurement are represented via the prob- 
ability vector Q — (q\,qi, - ■ ■), and P — (p\, pi, - ■ ■), and 
therefore 



Q = TP. 



(10) 



A general form is expressed as P„ = T„ „_iP„_i, where T ni „_i 
transfers vector P n -\ for state p„_i to P„ for state p n (n > 2). 
mjj could be called the quantum transition probability, and 
has the properties 



Y j m ji =hY u P 



iMfi = 1. 



(ID 



T is a lefty stochastic matrix, of which each column sums 
up to one and all elements are nonnegative. If the ranks of 
the two measurements are the same, T is square matrix. For 
Markovian case, by comparison with Chapman-Kolmogorov 
equation ([§}, the correspondence is qj ~ p(xi,t^), p-, ~ 
p(x 2 , t 2 \x u h), rriji ~ p(x 3 , t 3 \x 2 , h). 

Furthermore, the matrix T has the feature that its formula 
depends on the measurement steps and the quantum state. As 
we know the channel formula is fixed for a simple Marko- 
vian quantum dynamics, however, the matrix T varies gen- 
erally, which we treat as the manifestation of the property 



of contextuality. Open-system quantum dynamics is contex- 
tual generally. A context is formed by arbitrary two succes- 
sive measurements M n -\, M n , and the quantum state p„_2 on 
which the measurement M n -\ is implemented. Notably, we 
find that when the measurements in the chain reduce to pro- 
jective measurements, the matrix T will not depend on the 
system state, and thus the dynamics is "de-contexted", and 
then we say there is no contextuality anymore. For the con- 
struction above of the measurement-chain with two measure- 
ments, if they are projection-valued measurement (PVM), and 
binary for simplicity, namely, M = {\m + ){m + \, \m~){m~\} and 
N - {|n + )(n + |, \n~){n~\), the matrix T takes the form 



T 



[<m> + >| 2 \(m-\n + )\ 2 



\{m+\n-)\ 2 \(m-\n-)\ z 



(12) 



which is a doubly stochastic and symmetric matrix. If there 
exists evolution between the measurement processes, it will 
not be symmetric. It is interesting that the measurement trans- 
fer matrix formula relates to the scattering matrix (S-matrix) 
approach to particle physics, and the general form of T with 
POVM measurements could be treated as a generalization of 
the S-matrix. For the chain of PVM, if each measurement is 
the same, the matrix T becomes identity operator. 

The property of contextuality of quantum dynamics con- 
cerns the differences between PVM and POVM. In the chain 
of PVM, the information of the initial quantum state is lost, so 
in this sense, PVM is maximally destructive. On the contrary, 
POVM does not destroy the information of the initial quan- 
tum state completely at once, so there will be coherence and 
thus contextuality in the dynamics, until the system is even- 
tually decohered. We will analyze some examples in the next 
section to demonstrate the points. 



B. Examples: contextuality in quantum dynamics 

In this section we employ the quantum measurement-chain 
formula to analyze some standard physical models. First, 
let us consider the exponential decay process described as 
amplitude-damping channel. Starting from the state |^o) = 
a\g) + /3\e), the system decays to the ground state \g) for equi- 
librium. We model the decay process as a chain of POVM 
measurements, each of which takes time r, and the two Kraus 

operators are K x = ^ K 2 = |jj where y = 

e -r/T^ j j s ^ jjf e tmie Q f exc j]; ec [ statg | e ) [i] The system 

state after the «-th measurement is 



Pn = 



1 -p 2 y n a/3y n/2 
a/3y" /2 (3 2 y n 



(13) 



The measurement transfer matrix between the (« - l)-th and 
«-th measurements takes the form 



Trt,n— i 



( l-y"- 2 +y" 
<- +r - 

-r" 



l-y- 
y" 

\ l-y"~ 2 +y"~ 



1 



(14) 
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and the asymptotic form at equilibrium is Y a 







The 



ihji forms the weak-measurement transfer matrix T. It has the 
properties 



relations T^ n l = tr^^iP^^/tr^/Tip,,^) etc are 
used to calculate the elements of T„ „_i. It is clear that the 
formula of the measurement transfer matrix Y„„_i depends 
on the measurement steps. From its formula, we can see that 
the transition probability from K2 to K\ is 1, which intuitively 
manifests the roles of the two Kraus operators, namely, K2 
causes the decay of the excited state population, while K\ 
causes both decay and decoherence, and in the long-time limit 
the effect of K\ dominates. Although decay process also exists 
classically and some quantum decay processes can be mod- 
eled by classical dynamics, however, classical decay process 
is not contextual, since no coherence exists at least classically. 
As the result, we see that the quantum decay process is con- 
textual, and contextuality captures the primary properties of 
quantum dynamics. 

Next, we analyze the Stern-Gerlach experiment with se- 
quential measurements, where no contextuality exists. A 
beam of electrons will split into two branches after passing 
through the magnetic field. The effect of the magnetic field 
is a selective quantum measurement, i.e. the branches of the 
state are selectively distinguished via different spatial paths. 
Note the POVM in quantum decay process above is not se- 
lective, though. Suppose the sequential measurements are al- 
ternate along z-axis and jc-axis. Based on formula d!2t . the 
transfer matrix becomes 



T = (H 

\2 2 



(15) 



which is all the same for each step, and thus the dynamics is 
de-contexted. Still, there exists quantum property responsible 
for the Stem-Gerlach experiment, which is the superposition 
of quantum state. The state | f.i) will be a superposed state 
when expanded and measured in the basis of | 1 Z , J, z ). This 
example indicates that the notion of contextuality is different 
from superposition, yet both of them rely on the existence of 
quantum coherence. 



C. Weak-measurement transfer matrix 

In this section, we show that the formula of measurement 
transfer matrix can be generalized. Suppose the first mea- 
surement Ai on system po is weak, in the sense that we 
record some information of the system (observable) without 
changing the state significantly [31], the system state becomes 
po - Po, and then perform the second measurement N which 
provides the post-selection. The probability becomes 

qj = HNjNjpo) 
i 

= X tr(MM,p ) 

= ^jthjipi, (16) 



(17) 



T is a generalized lefty stochastic matrix, not all elements are 
real. It relates to weak value 1I3H1 . a generalization of eigen- 
value and expectation value, which even can take complex val- 
ues. Suppose the two measurements are performed to measure 
observable A. Then, in weak measurement, with probability 
qj we measure the weak value m J A which takes the form 



ti(N]NjAp ) 
lr(.V ; .V,p„) ' 



(18) 



with (A) = tr(Apo) = 2; to a'7>- The weak value m J A reduces 
to the original form for PVM case |3~lj]. Note that in equa- 
tion COD, n)N; can also be viewed as one observable with 

j J 

expectation value qj, then the weak transition probability ihji 
is just the weak value of NjNj pre-conditioned on the opera- 
tion m\M[. 

IV. DISCUSSION AND CONCLUSION 

In this work, the properties of open-system quantum dy- 
namics, namely positivity and contextuality, were studied. 
These two notions are defined from different standpoints (such 
as there are differences between nonlocality and entangle- 
ment). Quantum dynamics is distinctive from classical me- 
chanics and stochastic dynamics, but there is no unique con- 
cept that could capture the essence of quantumness or open- 
ness completely. From our study we can see that the property 
of contextuality is much more basic than complete positivity. 
We demonstrated the limitation of completely positive (CP) 
map. We verified that the CP map can be modeled physically 
by the unitary-universe dilation model, which provides an el- 
egant method to resolve the initial correlation problem. As 
a subset of physically positive map and hermitian map, CP 
map is not sufficient to describe all kinds of quantum dynam- 
ics. Therefore we develop the quantum measurement-chain 
formula, which is natural to describe the property of contex- 
tuality of quantum dynamics. The fine-grained structures of 
quantum operators (including density matrix, observable, and 
measurement etc) would result in realistic effects since quan- 
tum dynamics occurs not only in Hilbert space, but also in 
spacetime associated with quantum reference frame such that 
different basis and decompositions exist physically. Quantum 
contextuality directly relies on the fine-grained structure of 
quantum operator and quantum coherence. The traditional 
way to reveal contextuality is to employ inequalities of expec- 
tation values of non-commuting observables, which indeed re- 
lies on the fine-grained structure of operators too. The formula 
of measurement transfer matrix can be viewed as a new way 
to understand contextuality. 

The framework of using swap operation to resolve the ini- 
tial correlation problem could have further applications in 
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channel theory, and quantum control etc. Also, problems of 
quantum and semi-classical non-Markovian quantum dynam- 
ics with initial correlation could be studied by applying swap 
trick. The measurement, and weak-measurement transfer ma- 
trix formula can be used in measurement theory, and analysis 
of quantum-classical transition etc. 
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Appendix A: Heisenberg picture 

An affine map form associated with Bloch coherence vec- 
tor, which is geometric and more general than CP map, can 
describe any hermitian (not necessarily positive) dynamics. 
In generalized Bloch ball, the quantum state is represented as 
a coherence vector, whose evolution manifests the quantum 
dynamics. This geometric picture is mostly analogue of clas- 
sical mechanics, namely, the dynamics of classical position 
or velocity vector can also be expressed as affine map. The 
Hilbert space equipped with Hilbert-Schmidt operator inner 
product (A, B) = tr(A 1 B) is suitable to describe open-system 
quantum dynamics. For a qudit system, there exists a com- 
plete orthonormal operator basis {Fj, i — 0, 2, . . . , d 2 - 1 }, and 
tr(FjFj) — 6/j, with Fq — -^=1, the others traceless. The qudit 

j2_i 

system can be decomposed as p = (l+2 i=1 fiFd/d, with/; as 
the expectation value of observable Fj, The Bloch coherence 
vector is formed as a column vector F = (/i,/2, • • • ,fd 2 -i)- 
The open-system dynamics is then represented by the Bloch 
vector, i.e. the evolution of the d 2 - 1 parameters fi. Heisen- 
berg picture has one merit that there exists one well-defined 
form for any hermitian map 113 211 on a system observable 0$ , 
which is 

^{O s ) = Y l d a pKlO s Ki 3 , (Al) 

where d a p forms a hermitian matrix. Once we know the form 
of the map, the dynamics for each Fj then the expectation 
value fi can be derived. Another equivalent way is to rep- 
resent the map in the basis of Fj as a matrix %, with elements 
%j = tr(Fi®(Fj))/d. Note, the H matrix will relate to Choi's 
matrix C in one elegant way (ikffijl) = (ij\C\kl), via involu- 
tion. The channel matrix % takes the form IB3I1 




(A2) 



and the dynamics becomes affine map F' = RF + r, and p' = 
Hp. 

It seems it is better to employ the Bloch vector dynam- 
ics via hermitian map or affine map formulas in Heisenberg 
picture than the positive map formula in Schrodinger picture. 
The Bloch vector dynamics is the quantum analogue of classi- 
cal mechanics, however, the geometrical visualization is only 
proper for low dimensional case, e.g. qubit dynamics in NMR 
research. 

Below, we present the CP map form on system observable 
Os . For the u-u model, one environment is needed with the en- 
vironment observable 1 E . Note for the initial correlated case, 
the environment is bipartite (including the ancilla). Denote 
the global unitary as U. In order to find the OSR form for the 
system observable, we apply the fact that 

(O s ) = tr s (O s p , s ) = tr s (<y sPs ) 

= tx SE ((tfO s ® \ E U)p s ®p E ) 

= trsiJ^Me^Os 8 l £ [/)p s |e'XeV>) 

U 

= tr s (J] Me^Os 9 UU\e l )p s ), (A3) 



then the OSR form of the system observable 0' s is 



0' s =2^<e i l^ t O s ®li E l/|e , '> 

i 

^Yj^^'Os ®W){e i \U\e i ) 
ij 

-Y j K;p s K il . (A4) 

i.j 



with Kraus operator Kjj = yfXi{e-'\U\e'), consistent with forms 
d2]i. The above formula applies to both cases (with/without 
initial correlation) for the u-u model. A special form of her- 
mitian map can be realized, where the pure state vector \e l ) 
can be further decomposed as |e') = 2 ff c a \ai), the Kraus op- 
erator K/j becomes K/j = yfT, ^ aa > c a c* a ,(a'.\U\ai). Then it is 
direct to verify that when the environment observable reduces 
to a projector instead of identity, and the environment state 
becomes pure, the CP map reduces to the form of hermitian 
map ( lAlb . with separable coefficients d a p = c a c*„. 
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